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We propose a simple way to determine the periodicities of wave packets in quantum systems directly 
from the energy differences of the states involved. The resulting classical periods and revival times 
are more accurate than those obtained with the traditional expansion of the energies about the 
central quantum number n, especially when n is low. The latter type of wave packet motion occurs 
upon excitation of highly charged ions with short XUV or X-ray pulses. Moreover, we formulate 
the wave packet dynamics in such a form that it directly reveals the origin of phase shifts in the 
maxima of the autocorrelation function. This phenomenon, so far poorly understood since it is not 
prominent in the high n regime, becomes a dominant feature in low n wave packet dynamics. 

PACS numbers: 32.80.Aa, 32.80.Ec, 42.50.Md 



I. INTRODUCTION 

A wave packet (WP) comprising many energy eigen- 
states of a system spreads during its time evolution, only 
to reverse the spreading and reshape after a certain time 
called the revival time. Originally formulated in the con- 
text of highly excited electrons in atoms [l| following the 
first observation of wave packet collapse and revival 0, 
the phenomenon has been identified in a large variety 
of physical systems, such as in Gaussian WPs in quan- 
tum boxes [3| , the evolution of rovibrational nuclear WPs 
[1, HI (for a review see @), even at an attosecond time 
scale Q, revivals in a coherent photon field Q, or re- 
cently, the propagation of WPs in graphene under mag- 
netic fields 0. 

Quantum revivals have also been widely studied over 
the years in a more mathematical context, with emphasis 
on obtaining accurate analytical expressions [13, [Hi , for 
a review see [l^ . 

The multiple interferences of the wave packet com- 
ponents lead also to the so-called fractional revivals at 
divisors of the revival time, which have been predicted 
[l3l [l^ and observed in Rydberg wave packets [l^ or 
photon bouncing balls [l^. Fractional revivals can re- 
trieve information about the system even if it decays be- 
fore the first revival [TTj . They can be used for mapping 
the quantum phase of a molecular nuclear wave-packet 
in two-dimensional spectroscopy |18j , or even to factorize 
prime numbers [l^. It has also been shown that infor- 
mation entropy in position and momentum spaces can 
reveal the existence of fractional revivals [l^l • 

The study of quantum revivals in atoms has tradition- 
ally focused on the excitation of Rydberg WPs. The 
latter contain typically of the order of 10 highly excited 
states when they are centered on a principal quantum 
number 30 < n < 85 and excited by laser pulses with a 
duration of a few picoseconds [l|. Under these circum- 



stances the relative difi'erence between the energies of the 
states is small, and the standard approach of representing 
the energies of the contributing states by Taylor expan- 
sion about the central energy En provides good results. 

However, new FEL or high harmonic-based light 
sources are able to deliver short pulses of 100 atomic 
units (2.4 fs) duration and less in the VUV to X-ray 
regime [2l| . In particular the X-ray pulses will typically 
lead to a multi-electron wave packet of the valence shell of 
an atom by inner shell ionization. In order for the photo 
electron to form a bound wave packet while simultane- 
ously the core electrons remain in their ground state, such 
wave packets must be generated in highly charged ions, 
as produced, e.g., in state-of-the-art electron-cyclotron 
ion sources (ECR) [12 or EBITs [13]. The high ionic 
charge implies large electronic energy spacings so that 
only a few levels contribute to the wave packet, despite 
the large energy width due to the short pulse duration. 
Hence, the dynamics lays in between quantum beats of 
a few well defined energy levels [2j| and the traditional 
regime of Rydberg wave packets as discussed before. 

A good estimate of the dynamical range of a system, 
which goes from quantum beats to wave packet revivals, 
is the number of states in the WP, k. One can estimate 
this number by relating the energy spread of the laser 
pulse to the energy spread in the excited atom, which 
for a WP centered in state n is ~ kdEn/dn\n- An 
electric field of Gaussian shape f{t) = /oexp~^'"^* , 
where T is the FWHM of the pulse and /o its amplitude, 
will populate a wave packet with amplitude 



/o cxp 



-{E„-E-y-T^/ {8 log 2) 



(1) 



for each energy level En- Including all those states with 
probability > c implies an energy spread of AE = C/T, 
where C ~ 4(— In 2 Inc)^/^. For a hydrogenic ion of 
charge Z and spectrum En = — Z^/(2n^), equalling both 
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TABLE I: Number of states k contributing to a wave packet 
with up to an amplitude weight of c = 10"'^ when excited to a 
central state n by a pulse of duration T (FWHM). We assume 
an hydrogenic ion with energy spectrum En ~ —Z'^/{2in?). 
Parameters given are used in the figures of this paper as in- 
dicated. 



_k Z n T [fs] Figs. 

13 1 85 8000 1 

8 1 45 2000 4,5 

4 8 14 2 4,5 

4 12 7 0.1 6 

5 8 9 0.4 8 




Energy (a.u.) Energy (a.u.) 

FIG. 1: (Color online) Field intensities in the frequency 
regime, and energy levels for two regimes: (a) high 

n, long pulses, for the particular case oi Z = 1, n = 45, 
FWHM=2 ps, (b) low n, short pulses, for the case of Z = 8, 
n = 14, FWHM=2 fs. The red vertical lines correspond to 
n in each case. For the number of energy levels participating 
in the wave packet with an amplitude of more than 10~^, see 
table 1. 

values of AE gives 

k = Cn^/{Z'^T). (2) 

A few quantitative examples of the number of states k 
as a function of n and T for different atomic numbers is 
shown in Table 1 . Two cases of low and high n are shown 
in Fig. m 

Our goal is to treat the low n regime of electron ex- 
citation from the perspective of wave packet revivals. 
However, due to the large energy spacing the standard 
description of the WPs with a harmonic approximation 
about En becomes rather inaccurate. In the following 
we will show that one can replace the derivatives of the 
energies by exact differences of energy levels to obtain a 
description which is valid in the new regime of short, 
high frequency pulses. This approach is conceptually 
very simple and it is universally valid, also in the tra- 
ditional Rydberg regime. We can derive accurate values 
for the revival times in all situations. Moreover, the ori- 
gin of the shift in the position of the maxima around 
the revival time, previously noticed [l^, [25| but not well 
described, emerges naturally. 

The structure of the paper is as follows: in section HIl 
a form of the autocorrelation function is derived which 
allows one to explicitly read off the relevant periods. We 



also derive and explain why for the low excitation regime 
finite energy differences should be used to define these 
periods instead of energy derivatives. The periods, clas- 
sical, revival and fractional revival times, as well as a fea- 
ture which is important for low excitation wave packets, 
namely the time where the extrcma flip their sign, are 
explained in section Hill with examples from the high and 
intermediate excitation domain. Section IIVI points out 
the characteristic features of low excitation wave packets 
ending with a hydrogenic ion excited by a 100 attosec- 
ond X-ray (1.92 keV) pulse. Finally, for multi-electron 
atoms, the effect of the quantum defect is discussed inlVl 
and conclusions are drawn in IVII 

II. AUTOCORRELATION FUNCTION 

Without loss of generality but to be somewhat specific 
we consider a system which absorbs one photon from a 
laser pulse of length T in a transition from an initial state 
I a) to a final state |n) which, for simplicity, we assume to 
be an cigenstate of the system. Due to the finite length 
of the pulse, this process creates a wave packet which can 
be expressed in first order perturbation theory after the 
pulse as 

= i^|n)exp[-i£;„t]Af„„^^„, (3) 

n 

where Fn are the amplitudes defined in Eq. ([T]). The 
dipole matrix elements Man = ("■|'Z|5'q(0)) give the 
strength of photon coupling to the different eigenmodes. 
Atomic units are used throughout the paper unless oth- 
erwise stated. 

The modulus square of the autocorrelation function 

A{t) = |(a|*(t))|' = |c„nc™pexp[-z(S„ - E,n)t] 

n.in 

(4) 

with coefficients |c„p = |A'/q,„P|F„P is formally equiva- 
lent to PQ,_j./(t), the probability of going from a state la) 
to a state |/) in a pump-probe experiment involving two 
identical pulses with a time delay t. Therefore, Eq. @ 
can be measured [l| . 

One of the salient features of Eq. ^ is that this time- 
dependent function shows revivals and partial revivals of 
the initial wave packet. Clearly not all of the infinitely 
many terms in Eq. (|4]) do contribute equally, since the 
pulse excitation covers only a finite energy interval of the 
spectrum {En}, see Eq. ([T]). Hence, our goal is to rewrite 
the double sum such that a systematic truncation to the 
most important terms is facilitated, yet retaining in the 
most accurate manner the different periodicities which 
give rise to the revival phenomena. 

To reorder the contributions in the sum, we will make 
use of an expansion of the eigenenergies En around the 
mean excitation £Vr, which is well justified as long as the 
pulse excites a region of high density of states, so that 
the relative level spacing {dE / dn\n) / En ^ 1. 
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In a second step we will acknowledge the fact that for 
high frequency pulses the spectrum of highly charged ions 
is excited in a region of much lower density of states than 
in the traditional Rydberg regime. 



A. Regrouping the sum 

First we perform a Taylor expansion of the energy 
around the mean excitation energy £Vr, and split the ex- 
pansion into two terms with even/odd powers of n, so 
that the energy differences become (see appendix I A ip 

Sn = 6+ + (5,7 (5a) 



where 



n+n 



where 



and 



2tt 



.7 = 1 



27ry ■ 



1 1 
2^7! 



d^E„+n 



(5b) 



(5c) 



For the regrouping we assume that the dipole matrix 
elements Man do not vary much in the relevant terms of 
Eq. ((4]) and are therefore irrelevant for the summation. 
Moreover, the diagonal parts m = n will only contribute 
a non-oscillating background, so that the time-dependent 
part of Eq. (U]) simplifies to 



a{t) = ^ Wjnn CXp[-i{Sjn " Sn)t] , (6) 



with weights 



w„ 



(7) 



decreasing rapidly with increasing distance S\n\ from the 
central energy. Here a = V41og2/T is the pulse width in 
frequency domain, see Eq. ([1]). The idea is now to arrange 
the double sum Eq. © into groups of terms which have 
similar Gaussian weight Wnm, where we order the weight 
according to dm on \m\. From a basic element Wnm a- 
quadruple of elements is generated by the two symme- 
try operations and their concatenation: exchange of the 
indices (n, m) — > (m, n) (which amounts to complex con- 
jugation), and point inversion (n,m) — > {—n,^m). Note 
that m, n can be negative due to the shift of the indices 
by the value n. Hence we may rewrite Eq. © as 

X! WmnCOa[{Sm-Sn)t]+W-m-nCOs[{S^m-S-n)t 
0<n<m 

(8) 

Making use of the symmetry properties of the Sn from 
Eq. ((5a|) . we get (sec appendix IA2[) 

{Cmn C0S[((5+ - S+)t] COs[{Sm - S-)t] 

0<n<m 

+ Smn sm[{S+ - S+)t\ sin [((J- - 5-)t\} , 

(9) 



isiy + is-y + is+y + is-) 



(lOa) 



accounts for the even and odd exponents of the neigh- 
boring levels, and 



C„ 



J=0 



(2j)!(T 



4j 



j=0 



(2j + l)!a 



are the coefficients of the cosine and sine terms in Eq. © , 
respectively. 



B. Expansion to lowest order around the mean 
energy 

The final form of the autocorrelation function Eq. ([9|) is 
still an exact version of Eq. ([6]) . While of little practical 
use, it reveals the structure quite well, which contains 
products of two trigonometric functions. Their mean- 
ing becomes obvious if we truncate the Taylor expansion 
Eq. ([5a|) around En to lowest order. 



2-KrTn? 
To 



2nm 



(11) 



Keeping only products of derivatives to second order, i.e., 
in times Tj"^,r2~^ and Tj"^, reduces Eq. ([9]) to (see ap- 
pendix IA3|) 



a2{t) 



0<n<m 



-^(m +n ) 



(12) 



27r(m^ - n^) 



To 



27r(m — n) 



From this representation one can easily read off the stan- 
dard essential periodicities in the autocorrelation func- 
tion: the classical and revival periods. One can also di- 
rectly see how periodic maxima of a{t) turn into minima 
for certain time intervals, as we will show in the next sec- 
tion. However, as we will also see, a further simplification 
using the two energy values En±i next to the mean en- 
ergy En can be made to construct the autocorrelation 
function 02 (i). Notice that this is not an additional ap- 
proximation: the phase in the cosine in Eq. (jS)), which is 
exact, is {6m — 6n)t = (£W+m — E)n+n)t, SO the discrete 
energy differences are the ones which determine the pe- 
riodicity of the system. 



C. Expansion to lowest order using finite 
differences in the energy eigenstates 

The autocorrelation function in Eq. (jH]) contains energy 
derivatives (see Eq. ([5c)) ). While they arc usually taken as 
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continuous derivatives, we will consider finite differences 



El 



2 2 

{En+i - 2En + En-l) _ (<5i + 5^i) 



(13) 



The motivation for this conies from the fact that for the 
X-ray excitation regime, the few states which contribute 
to the autocorrelation function (see table 1) are not nec- 
essarily narrowly spaced in energy, so that a Taylor ex- 
pansion may not be very accurate. 

On the other hand, using finite differences implies that 
one includes the most important terms in the autocorre- 
lation function with exact phases. The most important 
pair is {m,n) = (1,0), since it gives rise to the terms 
with highest weights in the autocorrelation function, see 
Eq. ([7]). For this pair of indices, the phases in Eq. 
become E^/2 and EL (Eq. (IT3|)). Therefore, these in- 
dices contribute to the autocorrelation function with the 
phases 



2 cos 



Si - <5_j 



cos 



cos [{6i - 5Q)t] + cos [((5_i - (5o)i] , 



(14) 



which are the phases of the two lowest order terms in 
Eq. IH]) (note that Sq = 0). The amplitude is shghtly off 
due to truncating Cmn in Eq. ^ to the j = contribu- 
tion, see appendix lA 31 However, this does not affect the 
periodicities of a(t), which we will discuss next. 



D. Difference of periodicities defined by continuous 
derivatives and finite differences 

A natural question is how much the periods Ti defined 
by the continuous derivatives in Eq. ()5a|) differ from the 
periods Tj defined by finite differences. This can be best 
illustrated with an example, for which we take a hydro- 
genic Rydberg electron with spectrum E„ = —Z'^/{2ii?) 
excited to E^. The standard definition for the classical 
period, understood as the orbiting period of a classical 
electron, is derived from the energy derivative 



- _ - _ 27r _ 2TTn^ 



while from finite differences we obtain 



(15) 



47r 



2n(lf 



En+l — En-l 



TZ2 



(16) 



The relative difference between both expressions is 



Tr.~Tr. 



1 



9 4 5 



(17) 



which is noticeable for lower values of n: it is around 1% 
for n = lA (with a relative level spacing AEn/En — 
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FIG. 2: (Color online) Relative differences between the stan- 
dard approach (continuous derivatives) and our method, in 
logarithmic scale. Black line: classical period (Eq. (|17|l ). 
Green line: revival time (Eq. (|19|) V 



2/n = 0.3), but less than 0.1% for n = 45 (with 
AEn/En = 0.07). Similarly the revival times are 



Tr 
Tr 



To 



T2 = 



47rn'' 

"3Z2"' 

47rn^(n^ 



(18) 



1)^ 



Z2(3n2- 1) 



with a relative difference of 



Tr 



3n" 



1 



(19) 



This amounts to ~ 0.85% difference for rl = 14, and only 
0.08% for n = 45. Hence, both expressions merge for high 
quantum numbers, while for low quantum numbers with 
a larger spacing between energy levels a finite difference 
is a poor approximation to a derivative, see the sketch 
in Fig. [1] In our case, the energy differences are the 
"true" ones, so continuous derivatives are an approxima- 
tion, which fails for low quantum numbers, as illustrated 
in Fig. [2] For this reason, we will use in the following the 
Ti derived from finite energy differences instead of the Ti . 



III. FOUR ESSENTIAL PERIODICITIES OF 
THE AUTOCORRELATION FUNCTION 

A. Classical period and revival time 

The autocorrelation function Eq. ^ or its truncation 
Eq. reveal directly the classical period Tc = Ti, for 
whose multiples the second cosine factor in each term 
becomes maximal. These are the two fundamental time 
scales emerging from the first and second order energy 
differences of energy levels next to the centrally excited 
one, as shown above. 

A representative autocorrelation function for a radial 
Rydberg wave packet which exhibits these well known 
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FIG. 3: Autocorrelation function for (a) Z = 1, n — 45, 
FWHM=2 ps, and (b) Z = 8, n = 14, FWHM=2 fs. Full 
lines: versus Tc. Dashed lines: versus Tc. In (a) both are 
indistinguishable. The classical period Tc, revival time Tr, 
and the time for sign change Ts are shown with vertical lines. 



features is shown in Fig. |3l for the high- and low-n 
regimes. The classical and revival times are indicated. 
Note that results for periods Tj and Ti coincide to the 
accuracy of the figure in the high-n case (Fig. [Sji), while 
a discrepancy is visible in the Iow-tT case of Fig. [31d. 



B. The time of sign change 

At first glance, the wave packet in Fig. shows a re- 
vival at Tr/2, but a closer inspection reveals that this is 
not an exact copy of the wave packet, since its maxima 
are shifted by half a period (in between the vertical 
dashed lines). The revival at Tr is complete, with max- 
ima close to full periods of Tc- Although noticed before 
[ll,!!!], there has been no formulation of the autocorrela- 
tion function which would reveal the shift directly. This 
is achieved in Eq. ([5]) with the products of two trigono- 
metric functions and the cosine product dominating as 
revealed by Eq. The first cosine function with its 

slow period Tr can be interpreted as an envelope to the 
rapidly oscillating second one with period Tc- Hence, we 
identify a time for sign change Tg = Tr/4 at which the 
sign of the envelope changes, turning the maxima of the 
second cosine function into minima. This is the origin of 
the apparent shift in the position of the maxima. 



t(T) 



FIG. 4: Autocorrelation function for a Rydberg wave packet 
with n — 85 and a pulse of 8 ps (a). The different fractional 
revivals are shown with vertical lines. 



C. Partial revivals 

As it is well known, in addition to the revivals, there 
appear fractional revivals at times 



t 



pTr 



2(m2-n2) 



(20) 



where the phase of the first cosine in Eq. p2| is a mul- 
tiple of TT. However, as the values of n, m increase, the 
corresponding Gaussian weights decrease. Therefore, the 
higher the order of the fractional revival, the lower its 
intensity. An example is shown in Fig. |4l Following the 
philosophy of keeping only the largest contribution, one 
can directly see from Eq. (|12p how the partial revivals 
emerge by truncating the sum according to the dominant 
contributions of the Gaussian weights. 

From Table |TT] we expect that the largest contribu- 
tions will arise from the fundamental pair (n, m) = (1, 0), 
which leads to the fractional revival at Tr/2. The second 
largest contribution comes from the (n, m) = (2, 0) term. 



which has an exponent of n 



= 4. According to 



Eq. (pn)) this contribution (with p = 1) should give rise 
to a fractional revival at Tr/8. If the third largest term, 
the pair (n, to) = (2, 1) with an exponent of 5 is included, 
a new partial revival will appear as expected at T^/G, al- 
though with a very small weight. This evolution can be 
seen in Fig. [5] for Z 1, n 85 and FWHM = 8 ps. 
In Fig. [5^, the autocorrelation function emerges as more 
terms are added. The exact analytical value is shown in 
the top (thick black line). The periodicities due to the 
different terms of the sum are shown in Fig. [5}d. 

Overall, as can be seen from Eq. ([T^ , two factors deter- 
mine the maximum order of fractional revivals resolved 
in the autocorrelation function: the width of the pulse ax 
and the energy difference between the states adjacent to 
n, as well as the inverse (discrete) classical revival time Tc 
in the combination ax/Tc: The shorter the pulse (small 
or) and the smaller the energy differences (large Tc) the 
higher are tliese weights, and therefore more fractional 
revivals will appear in these cases. 

Higher values of n with many states involved will create 
higher orders of fractional revivals, as shown in Fig. Hl^a) 
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FIG. 5: (Color online) (a) Different truncations of the sum 
in Eq. for Z = 1, n = 85 and FWHM = 8 ps, including 

up to {n,m) = (1,0) (black), (2,0) (red), (2,1) (green), (3,0) 
(blue) and (3, 1) (orange). The exact analytical value A(t) 
(Eq. Q) is shown for comparison at the top of the upper 
figure (thick black). The expected position of the fractional 
revivals is shown with vertical lines (see text). In (&), the dif- 
ferent components are shown separately (without summation) 
for the same cases plus for the pair (3, 2) (in dark green). 



TABLE II: Exponents m'^ + n^ for the approximate weights in 
Eq. p2l) of contributions (m, n) giving rise to partial revivals 



Tr/{m^ ~n^), see Eq. (QOl 



pair 


exponent 


revival divisor 


(m, n) 


2 1 2 


(m^ - n^) 


(1,0) 


1 


1 


(2,0) 


4 


4 


(2,1) 


5 


3 


(3,0) 


9 


9 


(3,1) 


10 


8 


(3,2) 


13 


5 



for a Rydberg wave packet with roughly 13 participating 
states (see table 1) with Z = 1, n = 85 and FWHM=8 
ps. 



IV. WAVEPACKET DYNAMICS IN THE LOW 
EXCITATION DOMAIN 

We will now explicitly discuss the new regime of low 
quantum number wave packet dynamics. As sketched in 
Fig. [T}d, for the case of a highly charged ion, i.e. with 
an effective atomic number Z = % (Fig. [TJd), when the 
n = 14 level is hit with a pulse of width FWHM=2 fs, 
only fc ^ 4 levels are excited and the relative energy dif- 
ference between the levels is large, thus one would expect 
periodicities obtained from a Taylor expansion to be less 
accurate. 

The corresponding autocorrelation function (Eq. (jj])) is 
shown in Fig.[3]D, both against (full lines), and against 
Tc (dashed lines). The discrepancy is clear and becomes 
already seizable around the first revival, for which Tc de- 
scribes very poorly the periodicity and Tc remains almost 
exact. Also prominently visible arc the maxima located 
at half-integer multiples of Tc around Tr/2 and at integer 
multiples of Tc close to the revival time T^, as explained 
in Section Hirei 

The extreme case of a very low quantum number in the 
final state is shown in detail in Fig. [51 Here we consider 
a hydrogenic atom with an effective Z = 12 excited to 
n = 7, with an ultrashort pulse of 0.1 fs with only 4 states 
participating, a situation which one may also consider as 
a multi-quantum beat phenomenon. The X-ray photon 
energy required for a Is — >■ 7p transition in this system is 
ijj = 1.92 keV (we neglect the quantum defect). While the 
sign change of the extrema at Tr/4 is very clearly visible, 
revivals can no longer be identified. Also, the classical 
time Tc from a continuous energy derivative completely 
fails to predict the time scale (and therefore the position 
of the extrema) correctly after about 2 fs. 

In Fig.[7)D we see that the particular value of the dipole 
elements has a small influence on the autocorrela- 
tion spectrum in this regime. In this figure, the same 
autocorrelation function as in (a) is shown (black line) 
with all dipole elements Mi„ = 1, together with the re- 
alistic calculation for Min = i^jnf'l'^ (green line) [26j . 
The details of a(t) vary, but the position of the peaks 
remains overall the same. 



V. NON-RESONANT EXCITATIONS AND 
QUANTUM DEFECTS 

In a general excitation process, the photon energy may 
not be resonant with any particular electronic transition. 
In this case 



(21) 



We can identify this with an effective energy iJ^j+An, 
where < An < 1. As An grows, the central energy goes 
from En to £W+i , and the classical period evolves approx- 
imately with the third power of n -I- An (see Eq. P^ ). 
Notice that what really happens is that while the pairs of 
eigenvalues involved in the different beatings remain the 
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FIG. 6: (Color online) Autocorrelation function for Z = 12, 
n = 7, FWHM=0.1 fs. (a) Black line: a{t). Red line: 
cos(27rt/rc). Green line: cos{2-Kt/Tc). Blue line: cos(27rt/r^). 
The vertical lines are multiples (full lines) and half-multiples 
(dashed lines) of Tc. (b) Same function compared with the 
case of varying dipole matrix elements Min, see text (green 
line). 
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FIG. 7: (Color online) (a) Same as in Fig. [6l but for four 
photon energies, corresponding to effective values n — 7, 7.3, 
7.6 and 8. The values of Tc and 5Tc are drawn in each case 
with vertical lines, (b) Fourier transform of the cases in (a), 
using the two first terms in Eq. (|12[) . Abscissas are expressed 
in terms of time as f = 27r/a;. Vertical lines show the different 
periodicities, see text. 



same, the relative weights in Eq. ([TJ arc changing, and so 
is Tc. This is illustrated again for Z = 12 and FWHM = 
0.1 fs in Fig. [71 for four different photon energies: 1.919 
keV (Is 7p transition), 1.922 keV (n=7.3), 1.925 keV 
(71=7.6) and 1.929 keV {n=8). The autocorrelation func- 
tion and the Fourier transform are shown for the four 
cases in (a) and (b), respectively. In the latter case, only 
two terms in Eq. (jl2[) are used for simplicity, and ab- 
scissas are shown as a function of time. Vertical lines 
indicate the different periods, for Ta-b = 2Tr/{Ea — Eb). 
While n = 7 and 8 show only two relevant contributions, 
more terms appear in the spectrum at non-integer values 
of n. 

We have explored so far the hydrogenic case in which 
only one electron is bound to a nucleus of charge Z. How- 
ever, if the system contains more than one electron, we 
might consider a nucleus of charge Z + p and p core elec- 
trons. The wave packet is then subject to an effective po- 
tential with a Coulombic tail of charge Z, whose eigenen- 
ergies are given by 



Z^ 

2nl 



(22) 



where tIq, = n — a. The quantum defect a depends on 
the angular momentum, but is almost independent on 
the quantum number n for n 3> 1. The shift in the 
energies affects the periodicities of the system [23, [28j : 
With the period for the hydrogenic case (a = 0) and 
the multi-electron case (a 7^ 0) denoted by Tc and T", 
respectively, the relative shift due to the quantum defect 
a is given by 



Sn{a) 



1 



a 



Tc (772 _ 1)2 n„ 

1 -f ncji[2 — c? -\- 3a77 — 377^) 



77„(7T2-1)2 



(23) 



This shift accumulates in time and should assume the 
value 



^Srr(a) 



27i^ 



377^ 



1 



-2a 



(24) 



around T^. This prediction is valid even in the low 77 
regime, as is shown in Fig. [S] for Z = 8, 77 = 9 and a 
pulse length of 0.4 fs, using both, a = and 0.3. 
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FIG. 8; (Color online) Effect of the quantum defect on the 
autocorrelation function, for Z = 8, n = 9 and FWHM = 0.4 
fs, using a = (black) and a = 0.3 (green). Time is plotted 
in units of Tc, the classical period for the case with a = 0. 
The shift in the peaks around T,. is ~ —2a. 



VI. CONCLUSIONS 

We have explored the dynamics of electron wave pack- 
ets in the new regime of low atomic excitation which 
can be realized by short and intense laser pulses of high 
photon energy, e.g., from XFEL light sources. We have 
shown that the standard analysis of wave packet dynam- 
ics with characteristic periodicities carries over to the low 
excitation regime, but only when the periods are calcu- 
lated from differences of the two energy levels bracketing 
the main excitation level n instead of using continuous 
derivatives. Moreover, we have formulated the autocor- 
relation function in such a way that a new time scale 
appears as an important structural element, namely a 
quarter of the revival time Tg = Tj./A, where the sign 
of the extrema in the autocorrelation function flip sign. 
This is a prominent feature in low excitation wave pack- 
ets and explains the "shift" in the maxima noticed pre- 
viously in the high n domain. Hence, the analytical and 
numerical tools are now available to analyze experiments 
on wave packet dynamics with intense short light pulses 
of high frequency. 
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51164 during a visit of KITP at the University of Santa 
Barbara. PR acknowledges support in part from the 
MICINN projects FIS2010-15127 and JCI2009-03793, 
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Appendix A: Regrouping the sum in the 
autocorrelation function 



1. Terms in the Taylor expansion 

The spectrum can be expanded around £^1 



En+n ~ En = 5n ^ (5„ + (5„ 



(Al) 



The Taylor expansion is 
1 d^E„ 



E 



j ! dm^ 



(m - noY 



^ (2j - 1)! dm2J-i 



1 d'^^Er, 



^ (2j)! dm^j 
At m = n + no it is 



n—TiQ 

{m-no)^^ (A2) 



E (2j - 1 

i=i ^ 



1 



n+no 



)! dn2i-i 



^ 1 d^^E„ 
^(2j)! dn 



2j 



n=0 



(A3) 



so if we define 



then 



1 _ 1 d^En+no 

Yj ^ 2^ dn) 



n=0 



(A4) 



2-E 



I T20-I 

+ 



27r > ^ 



(A5) 



2. Final expression for a{t) 
An expansion of the cosines in Eq. ([8]) brings 



a{t)/2 



W^n cos[(<5+ -(5+)t]cos[(,5„-<5„)t] 



0<n<m 








- sinKCi 


-5+)t] sin[(*,„ - 








cos[((5i,„ - 5t„ 


)t] cos[((5_,„ 




-5+Jt]sin[(<5: 







Due to the symmetric properties of the energy differ- 
ences, it is = (5+ and 5Z„ = —5~, so 



The coefficients are 

where Wmn is defined in Eq. (jlOap . and 



(A7) 
(A8) 



(A9) 



(AlO) 
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Again due to symmetry reasons, w 

^ — TT). — r? ^ 



Wmn and 



g-«'™„/<T"g2«„„/<T" 



(All) 



and thus 
With this, Eq. (|A6|) becomes 

a(0/2 = ^ M/„„,[cos[(<5+ -<5+)t]cos[a- -<5-)t] 

0<n<m 

- - O^] sin[(5„ - S-)t] 



3. Lowest order of the expansion 



If we take the lowest order (j — 1) in the definition of 
(5+ and (5~ at Eq. (fTTj). then 



27r 
27r 



and 



+ sin[(e - S+)t] sm[{S-, - S-)t]\ . (A12) 
The terms with cosines gather in a term with coefficient 



(A18) 



= 47r^ 



-1 -^2 ^1 



2e-»™./<T' cosh 



2u„ 



(A13) 



while the terms with sines gather with a coefficient 



which at second order in the derivatives is 



An- 



(A20) 



26""''""/'"' sinh 



2w„ 



, . The lowest order of the sums in Cmn and Smn is j = 0, 
' so 



Therefore, 



(A15) 



Sri 



1, 

2w„ 



(A21) 



and 



2i+l 



(2j + l)!cr4j + 



#16) 



Since Smn = at second order in the derivatives, and 
using (Tt = 1/cr, we obtain Eq. (|12p . 



which are respectively 2Cmn and 2Smn, as defined in a2(t)/4^ 
Eq. (|10ap . This way we obtain Eq. o<n<m 



47r CT^ . 2 I 2\ 
+n ) 

e ^1 cos 



2Tr{m'^ ~ n^) 



a(t)/4 = J2 ^~ 

0<n<m 



C,nn C0S[(5+ - S+)t] C0S[((5„ - 5„ )i] 



27r(m — n) 



(A22) 



+5™„ sin[(J+ - S+)t] sm[{S-, - 6-)t] . (A17) 
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